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ABSTRACT 

This study deals with the spin wave spectrum in magnetic macrostructure (composed of 
two ferromagnetic materials) showing a 3D periodicity: spherical ferromagnetic grains 
disposed in the nodes of a simple cubic crystal lattice are embedded in a matrix with dif- 
ferent ferromagnetic properties. It is shown that the magnonic spectrum of this composite 
structure exhibits frequency regions forbidden for magnon propagation, and the energy 
gaps are found to be sensitive to the exchange contrast between the constituent materials 
as well as to the magnetization contrast. The widths of the respective magnonic gaps are 
studied as functions of parameters characterizing the magnetic structure. 
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INTRODUCTION 



Although the first study of electromagnetic wave propagation in periodic structures, writ- 
ten by Lord Rayleigh, was published already in 1887, it is only in recent years that 
photonic composites raised suddenly an extremely interest. The research in this field was 
initiated by the studies by Yablonovitch and John [TJE] published in 1987 and anticipating 
the existence of complete energy gaps in electromagnetic wave spectra in three-dimensional 
periodic composites, henceforth referred to as photonic crystals. These can be used for 
fabricating new optoelectronic devices in which the role of electrons, traditionally used as 
transport medium, would be played by photons |S1 IH E] • The so-called left-handed mate- 
rials (LHM), showing negative effective refractivity [HI [7], provide an example of periodic 
materials demonstrating how much the properties of this kind of structure can differ from 
those of homogeneous materials. Another type of periodic composites are structures com- 
posed of materials with different elastic properties; showing an energy gap in their elastic 
wave spectrum, such composites are referred to as phononic crystals jHl El ^] • Recently, 
attention has been focused on the search of photonic and phononic crystals in which both 
the position and the width of the energy gap could be controlled by external factors, such 
as applied voltage or magnetic field. Attempts have been made to create photonic crystals 
in which one of the component materials is a magnetic ^T] ^1 El El QE] • 
A magnetic periodic composite consists of at least two magnetic materials; the infor- 
mation carrier in such structures are spin waves. By analogy to photonic and phononic 
crystals, in which the role of information carrier is played by photons and phonons respec- 
tively, periodic magnetic composites are referred to as magnonic crystals. Studies of 2D 
magnonic crystals have been reported [13 UHl CHI EUl 122 ESI , with scattering centers in the 
form of "infinitely" long cylinders disposed in square lattice nodes (cylinder and matrix 
materials being two different ferromagnetics), and the anticipated gaps were found indeed 
in the respective spin wave spectra. Further research was focused on magnetic multilayer 
systems, which can be regarded as ID magnonic crystals P I 1 ^ 1 ^ l?7| 1 ^ 1 ^ IHHl I!TT 1 1!?^ . 
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In this paper, we present numerically calculated magnonic band structures of three- 
dimensional magnonic crystals. Due to the complexity of the problem, only the simplest 
model of 3D magnonic crystal is considered here, represented by a system of ferromagnetic 
spheres (which act as scattering centers) disposed in the nodes of a simple cubic crystal 
lattice and embedded in a different ferromagnetic material (matrix). Both the exchange 
and dipolar interactions are taken into account in our calculations based on the plane 
wave method and using the linear approximation. 

THEORY OF 3D MAGNONIC BAND STRUCTURE 

Let's consider an ideal periodic structure consisting of spheres of ferromagnetic A and 
embedded in a matrix of ferromagnetic B. The spheres are assumed to form a 3D periodical 
lattice of sc type (Fig. [T^). A static magnetic field, Hq, is applied to the composite along 
the z axis, and assumed to be strong enough to saturate the magnetization of both 
materials. The lattice constant is denoted by a; the filling fraction, / = 4/37TR^a~^ , is 
defined as the volume proportion of material A in a unit cell. 

a) b) 



A 




Figure 1: a) The 3D periodic structure studied in this paper; the structure consists of 
ferromagnetic spheres of material A embedded in a matrix of material B (materials A 
and B have different magnetic properties); the spheres are disposed in the nodes of a sc 
lattice, (b) The Brillouin zone corresponding to the considered structure, and the path 
MVXMR (highlighted) along which the magnonic band spectra are calculated. 



Ferromagnetics A and B are characterized by two material parameters: the spontaneous 
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magnetization {Ms^a and Ms^b), and the exchange constant [Aa and Ab)', both these 
parameters depend on the position vector r = {x,y,z): 

Ms{f) = Ms,B + {Ms,A-Ms,B)S{f^, 
A{r) = AB + {AA-AB)Sir}, (1) 

the value of function S{r) being 1 inside a sphere, and beyond. 

In the classical approximation spin waves are described by the Landau-Lifshitz (LL) 
equation, taking the following form in the case of magnetic composites (with damping 
neglected): 

dM{f,t) 



dt 



7/ioM(r,t) X H,ff{r,t), (2) 



where magnetization M{f, t) is a function of position vector r and time t; Heff{f, t) stands 
for the effective magnetic field fHlCSl 12011^ acting on magnetization M{f,t): 



Heff{r,t) = Hoz + h{f,t) + — ( V ■ -^V ) M{f,t); 



(3) 



z is the unit vector along the z axis; /i(r, t) is the dynamic magnetic field due to dipolar 
interactions; the third component represents the exchange field. The magnetization vector 
can be represented as the sum of its static component, M^z, which is parallel to the applied 
field, and its dynamic component, m{f,t), lying in plane {x,y): 

M(f, t) = Mgz + m{f, t). (4) 

The dynamic dipolar field, h, must satisfy the magnetostatic Maxwell equations: 

V X h{r) = 0, 

V ■ (hir) + m(f)) = 0. (5) 

In a magnonic crystal, the position-dependent coefficients in (jHl), i.e. Ms and A, are 
periodic functions of the position vector, which allows us to use, in the procedure of 
solving the LL equation defined in (j21), the plane wave method, described in detail in our 



earlier papers IT^ (dealing with two-dimensional magnonic crystals). Following this 
scheme, we proceed to Fourier-expanding all the periodic functions of the position vector, 
i.e. the spontaneous magnetization, Ms, and parameter Q defined as follows: 
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Q = -TTTFT- (6) 



The dynamic components of the magnetization can be expressed as the product of the 
periodic envelope function and the Bloch factor exp{ikf) {k denoting a 3D wave vector); 
the envelope function can be transformed into the reciprocal space as well. Including all 
the expansions into (j21) and (jSJ leads to the following infinite system of linear equations 
for Fourier coefficients of the dynamic magnetization components, 'fn^^{G) and myj:{G): 

iflm^kiG) = niyj^iG) + ^ ^ Ms{G - G ) 

Qi -no|fc + iJ" I 

+ EE[(^+ ^') • {k + G") -{G- G") ■ (G - G')]Ms{G - G")Q{G" - G')m^g(G'), 

G' G" 

- E E[(^ + G') ■ {k + G") - (G - G") ■ {G - G')]Ms{G - G")Q{G" - G')m^j:{G'), 

G' G" 

kx, ky and Gx, Gy denoting the Cartesian components of the wave vector k and a reciprocal 
lattice vector G, respectively; a new quantity introduced in ((Tj) is fl, henceforth referred 

to as reduced frequency: 



The Fourier coefficients of spontaneous magnetization Ms and parameter Q are calculated 
from the inverse Fourier transformation; in the case of spheres the resulting formulae read 
for Ms as follows: 

f Ms,^/ + Ms,b(1-/), forG = 

and similarly for Q; R is the sphere radius (Fig. ^). 



Obviously, the numerical calculations performed on the basis of (|7j) involve a finite number 
of reciprocal lattice vectors G in the Fourier expansions; however, we have made sure that 



the number used is large enough to ensure good convergence of the numerical results. As 
indicated by an analysis performed, a satisfactory convergence is obtained already with 
343 reciprocal lattice vectors used. 




Figure 2: The magnonic band structure found numerically for the sc lattice-based com- 
posite (iron spheres embedded in a YIG matrix). The spin wave energy branches have 
been calculated along the path (in the first Brillouin zone) shown in Fig. ^ for the filling 
fraction /=0.2. The lattice constant value is assumed to be a = lOOA, and the applied 
magnetic field value is HoHq = O.IT. 



NUMERICAL RESULTS 

The 3D magnonic crystal studied is a magnetic composite consisting of ferromagnetic 
spheres (of material A) disposed in the nodes of a sc lattice and embedded in a different 
magnetic material (B), referred to as matrix (Fig. Q). The corresponding magnonic 
band structure will be calculated along path M = 7r/a(l, 1,0) — > F = 7r/a(0,0,0) 
X = 7r/a(0, 1,0) — i> M = 7r/a(l,l,0) = 7r/a(l, 1, 1) in the nonreducible part 

of the first Brillouin zone (see Fig. Q). Iron (Fe) and yttrium iron garnet (YIG) are 
chosen as component materials A (spheres) and B (matrix), respectively, in the studied 
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example. As established in our earlier studies P^ IT^ (in the case of 2D magnonic crystals) 
such composition, involving a substantial contrast of the magnetic parameters between 
the component materials (Ms^Pe = 1-752 10^Am~\ Ms^ig = 0.194 10^ Am~\ Ape = 
2.1 10"^^Jm~^, Ayig = 0.4 lO^^^Jm"^), leads to opening large energy gaps in the spin 
wave spectrum. Also in the spin wave spectrum obtained for the sc lattice-based 3D 
magnonic crystal considered here, two wide energy gaps are present, one between bands 
1 and 2, the other between bands 5 and 6, as shown in Fig. |21 The existence of the gap 
signifies that spin waves having frequency values within the gap cannot propagate through 
the composite. The following parameter values have been assumed in our calculations: 
crystal lattice constant a=100A, applied static magnetic field ^qHq = O.IT, and filling 
fraction /=0.2. 




0.1 0.2 0.3 0.4 0.5 60 80 100 120 140 

Filling fraction,/ Lattice constant, a 



Figure 3: (a) The magnonic energy branches in the Fe/YIG simple cubic magnonic crystal 
plotted versus the filling fraction. The lattice constant value is assumed to be a = lOOA, 
and the applied magnetic field value is /^o-f^o = O.IT. (b) Magnonic branches plotted versus 
the lattice constant a of the same sc structure; the filling fraction is fixed at / = 0.2. A 
frequency gap (the white region) moves down with increasing a. 

Let's examine the effect of the filling fraction on the width of the spectral gaps. Figure 
Eh shows the magnonic bands plotted against the filling fraction in the sc lattice-based 
structure (Fe spheres embedded in YIG). The plot was obtained through projecting the 
magnonic band structure calculated along MTXMR path (for a fixed filling fraction 
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value) on the reduced frequency (Q) axis; afterwards the procedure of projection was 
performed repeatedly for consecutive filling fraction values. The maximum width of the 
first gap is found to occur at /=0.27, reaching value Ai7=63.02. Unlike in ID and 2D 
magnonic crystals fH^ IT?H , no oscillatory variations of the energy gaps with the filling 
fraction are observed in the considered 3D composite. Figure Eb shows the computed 
magnonic band structure plotted against the lattice constant a for / = 0.2. The gap 
center (Qq) is found to descend, and the gap itself AQ (the white region) to narrow as 
the lattice constant increases, in such a way that the reduced gap width AQ/Qq remains 
almost constant. 

Let's proceed now to the role played by magnetic parameters characterizing component 
materials. Figure 0^ show spin wave energy branches plotted against the contrast between 
the spontaneous magnetization values in materials A (spheres) and B (matrix); this 
contrast, defined as ratio Ms^a/Ms^b, will be henceforth referred to as magnetization 
contrast. The computations have been performed for a fictitious matrix material whose 
spontaneous magnetization and exchange constant values, Ms^b and Ab, respectively, are 
fixed at values close to those in YIG; the exchange constant value in the sphere material 
is assumed to be Aa = 2.1 10^^^ Jm~^ . The other parameters are fixed as well: the 
lattice constant value is a=100A, the filling fraction /=0.2, and the applied magnetic field 
fioHo = O.IT. It is seen that the gap width is maximal for the magnetization contrast value 
13, which is not far from the magnetization contrast between iron and YIG {Ms^Fe/ ^syiG 
= 10.82); note also that the required value of magnetization contrast necessary for the 
energy gap to open is greater than 2. 

We shall now examine the role of the contrast between the exchange constant values in the 
component materials. The results of the respective computations are depicted in Fig. 
showing magnonic branches versus the exchange contrast, defined as the ratio {Aa/Ab) 
of exchange constant values in materials A (spheres) and B (matrix); the only variable 
in the computations was the exchange constant in material A. All other structural and 
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Figure 4: (a) Magnonic spectra plotted versus the magnetization contrast for a sc 
lattice-based structure. The exchange constant value in the spheres is assumed to be 
= 2.1 10~^^Jm~^. (b) Magnonic branches plotted versus the exchange contrast; 
the spontaneous magnetization in the spheres is assumed to be Ms^a = 2.1 10^ Am~^. 
The assumed spontaneous magnetization and exchange constant values of the matrix are: 
Ms,B = 0.194 lOMm-i and Ab = 0.4 lO^^Vm-^; the lattice constant a = lOOA, the 
filling fraction / = 0.2, and the applied field HoHq = O.IT. 



material parameters are fixed at the following values: Ms^a = 1-752 10^ Am ^, Ms^b = 
0.194 10^Am-\ Ab = 0.4 10"" Jm-\ a = lOOA, fioHo = O.IT. From the Fig. Ho one can 
infer that the exchange contrast is not an indispensable factor for the opening of complete 
energy gaps in 3D magnonic crystals. 

CONCLUSIONS 

We demonstrated that the 3D magnonic crystal with a sc lattice-based structure exhibits 
complete energy gaps, whose appearance is due to two magnetic contrasts: (1) the ex- 
change contrast and (2) the magnetization contrast. Gaps are found to exist only when 
the spontaneous magnetization contrast is greater than 2. However, the gap width can 
be also controlled through adjusting the respective structural factors: the filling fraction 
value and the composite lattice constant. 
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